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Abstract 

We discuss the neutrino mass matrix which predicts zero or small values of [ V13 j in the 
MSSM and found the inequality, sin 2 2# 12 < sin 2 29 & , where sin 2 2# 12 is the mixing angle 
at Mr scale and sin 2 26® is the value determined by the solar neutrino oscillation. This 
constraint says that the model which predicts a larger value of tan 2 # 12 at Mr than the 
experimental value is excluded. In particular, the bi-maximal mixing scheme at Mr scale 
is excluded, from the experimental value tan 2 6 & < 1. In this model, | V13 1 and a Dirac 
phase at mz which are induced radiatively may not be small. 
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1. INTRODUCTION 



The SuperKamiokande data has discovered the neutrino mixing between and v T 
from the atmospheric data|I[. Now the SNO data[0 together with SuperKamiokande 
data|| have solved the solar neutrino puzzle and pinpointed a solution among four so- 
lutions. That is, its origin is mainly due to the mixing between the v e and and the 
LMA solution is the most probable one. They are summarized as 



tan 2 # ~ 0.34 , Ami ^ 5 x 10~ 5 eV 2 



sin 2 29 



atm 



A ™atm - 



3 x 10" 3 eV 



(1) 



where 8 & and # atm are mixing angles which appear in the solar and atmospheric neutrino 
oscillations, and in effect are mixing angles between the 1st and the 2nd and between the 
2nd and the 3rd mass eigenstates, respectively. Am Q and Am 2 tm are the squared mass 



2 

atm 



m 2 — m\ | , where m, is the mass of 



differences defined by Ami = m\ ~ m \ an d Am 
the i-th mass eigenstate of neutrinos. Usually, the sign convention, Am^ = m\ — m\ > 
is taken, in which case the result from the SNO-SuperKamiolande data|2|, [3[ favors the 
normal side, tan 2 # < 1, and disfavors the dark side, tan 2 ^© > 1 @. Another important 
information from the CHOOZ data|| gives a severe upper limit for [ Vxs| 



V 13 \ < 0.16 



(2) 



where V\3 is the element of the MNS neutrino mixing matrix representing the mixing 
between the 1st and the 3rd mass eigenstates. If we combine these information, the 
neutrino mixing matrix is approximately written by 



V 



Sr. 











'atm 



3 '«M 



(3) 
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where we included two Majorana CP violation phases^, ^|, ^ in the mixing matrix which 
play an important role for the neutrinoless double beta decay 0. 

Among the above experimental information, the most mysterious point is the puzzle 
why |Vi3| is so small in comparison with other mixing angles, 9 Q and # a tm- If it is really 



1 



small, we have to find out the reason for it. For a small quantity at the low energy scale, 
the naturalness usually asserts that it is zero at the higher energy scale, because it is 
quite hard to reproduce such a small quantity at the low energy scale. 

In this paper, we consider a possibility that | V13 1 = at the energy scale where the left- 
handed neutrino mass is induced by the see-saw mechanism. There are many advantages 
to consider this possibility. (1) The small value of | V13I is naturally explained because it is 
induced by the radiative correction. (2) This scenario may be realized in some theoretical 
models at Mr scalefllO], [11], [12], [13|]. (3) The Dirac CP violation phase is induced by the 
radiative correction. 

Now, we consider the neutrino mass matrix which predicts | V^3 1 = at Mr scale, 
in the diagonal basis of the charged lepton mass matrix. This mass matrix contains 
only seven parameters, three neutrino masses, two mixing angles and two Majorana CP 
violation phases, and thus there is no Dirac CP violation phase at Mr scale. This may 
gives an possibility that the Dirac CP violation phase which appears in the neutrino 
oscillation may be related to Majorana CP violation phases, which may be related in 
the leptogenesis, since in our model, two Majorana phases are associated with phases 
of neutrino masses and they may well have some relation with phases from the heavy 
right-handed Majorana mass matrix. 

This paper is organized as follows: In Section 2, we explain our model and the frame- 
work of neutrino mass matrix. The radiative correction is taken into account and the 
mass matrix is diagonalized analytically to connect the parameters at Mr scale and the 
present experimental scale, m^. In Section 3, the general feature of our result is explained 
and the predictions are given. In Section 4, by using analytic result, numerical analysis is 
made on the induced size of [ X^x3 [ , the Dirac CP violation phase, 5, and the effective mass 
of the neutrinoless double beta decay. The discussion on the absolute size of neutrino 
mass is given. Summary and discussions are given in Section 5. 
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2. THE MODEL 



We consider a class of left-handed neutrino mass matrices which gives V13 = 0, where 
V is the neutrino mixing matrix. We assume that this mass matrix is derived by the 
see-saw mechanism according to the SUSY GUT scenario at the right-handed neutrino 
mass, Mr scale and evolves following the renormalization equation for MSSM to the Z 
boson mass scale, mz- In this model, V13 at mz scale is induced by radiative correction. 
We examine the size of V13, the Dirac CP violation phase, 5, Majorana CP violation 
phases, 01 m and Pm, and the effective mass of the neutrinoless double beta decay, (m v ). 
(a) The mass matrix at Mr 

The mass matrix which gives V13 = is generally expressed in the diagonal basis of 
the charged lepton mass matrix as 



where 



m u {M R ) = OD v O ri 



D v = diag(M 1 ,M 2 e iao ,M 3 e i/3 °) 



(4) 



(5) 



with Majorana phases, «o and (3q. O is the mixing matrix at Mr scale and is given by 

(6) 



(\ » I" n ^ 



o 



C23 — s 23 
s 2 3 c 23 



C12 -S12 
S12 C12 
1 



with Cij = cos 8ij and = sin^. In the following, we use only as an angle at Mr 
scale. 

(b) The neutrino mass matrix at mz 

In MSSM, the neutrino mass matrix at mz is given byfl4[ 



mvijnz) — diag(l, 1, a)OD u O diag(l, 1, a 



(7) 



where a is defined by 



-^(l+tan 2 /3)(m T /v) 2 



a = l-e = 



M, 



< 1 . 



R 



(8) 



Here m T is the r lepton mass, v 2 = t> 2 + v\ and tan/3 = v u jv& with V{ being the vacuum 
expectation value of MSSM Higgs doublet (Hj)(i = u,d). 

In order to estimate e, we assume the right-handed mass scale, M R and the region of 

tan (3 as 

M R = 10 13 (GeV) , 2 < tan/3 < 50 . (9) 

Then, with m z = 91.187(GeV), m T = 1.777(GeV) and v = 245.4(GeV), we find 

8 x 10~ 4 < e < 5 x 10~ 2 . (10) 

(c) Masses and the mixing matrix at mz 

The effect of the radiative correction to neutrino mass matrix has been discussed 



by many authorsJT^, 15], |T6j and the following is known. (1) The mixing angles are 
stable for the case of the hierarchical or the inverted-hierarchical neutrino mass scheme, 
rrii ^ m 2 ^ m 3 or m 3 mi m 2 . (2) The instability occurs for mi ~ m 2 . (3) The 



Majorana phases in neutrino masses may play an important rolc[|T6 . 

Since the stable case is well analyzed, we focus on the unstable case. That is, we 
consider the following mass relation holds at Mr scale, 



Mi ~ M 2 , < A 2 i < [A31I 
< A 2 i < Ml , eMl ~eMj < |A 3 i| , (11) 

where A 2 i = M| — M 2 , A 31 = M| — Mf and we chose the convention, A 2 i > 0. The 
diagonalization of the neutrino mass matrix is made analytically and the derivation is 
given in Appendix. 

In the following, we summarize the result derived in Appendix. As for neutrino masses 
themselves, corrections are small and of order eM,, because we are considering the sit- 
uation where M 2 ~ M| ^> A 21 ~ eMf . Thus, neutrino masses at Mr and mz can be 
considered to be the same. 

Tn% — Mi . (12) 
The radiative correction gives effect to the mass difference between mi and m 2 , 



while the mass difference between the 2nd and the 3rd receives only a negligible effect, 
so that 



A 2 2 2 2 2 a 

Am atm = m 3 - m 2 - m 3 ~ m l - A 31 • 



(14) 



In the above, we required Am^ > in accordance with the common experimental analysis 
which gives tan 2 9 Q = 0.34 < 1. As for mixing angles, the radiative correction does not 
give any effect to the mixing between the 2nd and the 3rd mass eigenstates either. That 
is, 

0atm = 23 . (15) 

Thus, in the following, we use m ; for Mj except for the discussion of the mass difference 
and # atm for 6*23 in order to express 9 Q , | V13 1 , 5, cum and (3m in terms of observables at 
the low energy as possible as we can. 

The MNS mixing matrix which is given in Eq. (A. 15) is expressed as 



V 



V 

where atm = 9 23 , 9 Q is 



Cq —Sq 
s c atm c c atm 
S©Satm Cq S a tm 



IT/ \c~ M>\ ( 

I M3 e 



Satm 
Catm 
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ta M 



(16) 



S12C + ci 2 se 



ia /2\ 



Cq = \C 12 C- S 12 Se iao/2 \ , 



with 6 defined by 



sin 29 



4em 2 sf tm sin 29 12 cos ^ 



The induced mixing element, |Vi 3 | is given by 



V 13 \ = 



emim 3 sin 29\ 2 sin 26 ) atm sin 



Am 2 



(17) 



(18) 



(19) 



atm 



A Dirac CP violation phase, 5, and two Majorana CP violating phases, au and (3m are 

s = 6+6-f + f -A, 

do 



olm — £2 — £1 — ^~ ' 







M 



6- 



00 

2 ' 



(20) 



with 

fi = arg(ci 2 c - s 12 se- iao/2 ) , 

6 = arg( Sl2 c + c 12S e ia °/ 2 ) . (21) 

In the mixing matrix, # 12 and two Majorana phases, cto and /3q are only parameters 
defined at Mr. All other parameters are expressed by physical quantities at mz- 



3. GENERAL FEATURES 

As we explained, we take the convention AjBq > for which the result from the 
SNO-SuperKamiokande data requires that the mixing angle should be in the normal 
side, tan 2 # = 0.34 < 1@, §. Also we take the convention, A 21 = Mf - M? > 0. 

(a) The solar mixing angle 

Here, we discuss the relation between 6*12 defined at Mr scale and 9® defined at mz, 
the value from the solar neutrino oscillation data. We parametrize the solar neutrino 
mixing angle as 

tan 2 9 & = jj-^p , or sin 2 29® = 1 - p 2 , (22) 
and then p > is required to guarantee tan 2 9 & < 1. From Eq. (17), p is given by 

CtQ 

p = cos 29\2 cos 29 — sin 2# 12 sin 29 cos — 

= cos 2# 12 cos 29 - 2h sin 2 2# 12 cos 2 y , (23) 

where we used sin26 l defined in Eq. (18) to derive the second line and we defined the 
positive quantity h to avoid the complexity of equation, 

»= 2j £t- < 24 > 

From p > together with h > 0, we find cos26' 12 cos 29 > 0. Now we look carefully the 
equation for Am^ > in Eq. (13). With A 21 > together with the above condition, 
only consistent choice is 

cos2fl 12 >0, cos2fl>0. (25) 



Now that the sign of cos 29 is fixed to be positive, we can eliminate cos 29 in Eq. (23). 
Thus, we can express tan 2 9 Q in terms of 9± 2 , «o and h, 

|~ tan 2 ^ = cos 29 12 \ ll-4h 2 sm 2 29 l2 cos 2 ^ - 2h sin 2 29 12 cos 2 ^ . (26) 
1 + tan 9q V 2 2 

This is the equation which relates the angle 9\ 2 at Mr scale and the solar mixing angle, 

Next, we solve the above equation with respect to cosa an d find 
1 

2b 12 

By requiring < 1 + cos a < 2, we obtain 



1 + cos a = \ (- 1 cos 29 Q \-h cos 2 29 12 + cos 29 12 Jh 2 cos 2 29 l2 + 2| cos29 Q h + l) (27) 

ft, sin 2fto V * / 



—5 - S[n ' 26& 9 . 9/A ^ < sin 2 29 12 < sin 2 2# , (28) 

sin 2 2^ + (|cos2^ |+4e S 2 tm (m 2 /Am|)) 2 " " °' 1 ' 

where we used the expression of h in Eq. (24). This inequality shows the region of sin 2 29\ 2 
at Mr where the experimental value sin 2 2# at mz is realized. 

Before discussing the meaning of these equations, it should be mentioned that the 
above result is valid as far as |Vi 3 (M R )| <^ \ sin0 12 | is satisfied, where sin# 12 and |Vi 3 (M R )| 
are quantities at Mr scale. 

Now we discuss the meaning of the inequality. Firstly, we comment that the equality 
sin 2 2# 12 = sin 2 26' (stable case) holds in several cases, (1) the hierarchical mass case, 
mi m 2 , where m\ ~ AtUq, (2) the small e case which corresponds to the small tan/3, 
(3) the special case, «o — tt, even with m\ ~ m\ 3> Amg. 

One of the most important observation will be the following: The inequality sin 2 2# 12 < 
sin 2 26 l holds for most of models, because physically feasible models have to predict 
small values of |Vi 3 |. Thus, we can generally say that the model which constructed at 
higher energy scale such as Mr must have tan 2 9\ 2 less than or equal to the experimental 
value, tan 2 9 & = 0.34. The bi-maximal mixing scheme which is realized at Mr is not 
acceptable from the present experimental data. This may give a big obstacle for model 
building, because the model should predicts the experimental angle which does not have 
any particular meaning in the stable angle case. On the other hand, for the unstable 
case, the model needs to predict smaller value at Mr and the radiative correction lifts 



the value to the experimental one, by the interplay among neutrino mass, tan j3 and the 
CP violation angle, a - 

(b) The size of the induced | V13 1 

The induced |Vi 3 | is given in Eq. (19). Since it is proportional to e(m 1 m 3 /Am^ tm ), its 
value is suppressed by Am|/Am 2 tm , in comparison with corrections to the mass squared 
difference for the solar neutrino mixing and the solar neutrino angle. If m 3 /mi > 1, some 
enhancement is expected. 

(c) The CP violation angles 

The Dirac CP violation phase, 5 is induced from two Majorana phases. Since a is 
deeply involved in determining the solar mixing angle, 5 aside from (3 can be determined. 
Thus, we define 

<5i = 5 + A) = 6 + 6-| + y , (29) 

for which we analyze numerically in the next section. We hope that the knowledge of the 
phase (3 may be derived from the information from the leptogenesis. 

(d) Neutrinoless double beta decay 

With mi — m 2, the effective neutrino mass for the neutrinoless double beta decay in 
this mode is simply given by 



(m u ) = 



~ mi I (c 12 c - s 12 se~ iao/2 ) 2 + (s 12 c + c 12 se iao/2 ) 2 e- ia ° \ 



= mi A / 1 - sin 2 29 12 sin 2 y , (30) 



where we neglect rri3V 2 3 , because |Vi 3 | is small. 



4. NUMERICAL ANALYSIS 



Physical quantities, tan 2 # , AttIq, sin 2 26> atm , Am 2 tm , |Vi 3 | and (m„) are invariant 
under the change of «o to — ao- The quantity 5\ changes to —Si under the exchange of 
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a to — a . In the following, we confine the region of a to be < a < tt to discuss 
above quantities numerically. 

The radiative correction is proportional to e, which is a rapidly increasing function of 
tan/3 as seen in Eq. (8). Therefore, the effect is smaller for smaller value of tan/5. In the 
following, we consider two cases, tan (3 = 50 and tan f3 = 20. For the numerical analysis, 
we use the experimental data given in Eq. (1). 

(a) The angle sin 2 29 12 

As we see from Eq. (26), the solar angle is determined by a Majorana phase a , the 
neutrino mass m\ and the mixing angle sin 2 2# 12 at Mr scale. Therefore, when we give 
the value of sin 2 29 Q , three parameters are constrained and the contour curve for a given 
«o is drawn. In Fig. 1, the contour plot of a in the sin 2 2#i 2 and mi plain is shown for 
tan/5 = 50 with tan 2 9 G = 0.34. The wide values of sin 2 29 u are allowed which may be 
seen from Eq. (28). 

A particular feature is that the most of region corresponds to n/2 < a < n. That is, 
if we choose the value of oiq in this region, almost any value of sin 2 2# 12 at Mr scale can 
reproduce the experimental solar angle with an appropriate choice of mi, which should 
be greater than, say, 0.02 eV. If the mass mi < 0.01 eV, the sin 2 26* 12 is stable and should 
reproduce the solar angle precisely at M R scale. 

(b) The induced value of |Vi 3 | 

As we see from Eq. (19), | j depends on four parameters, cto, mi, m 3 and sin 2 2#i 2 . 
Therefore, we define | V13 1 (mi / 777.3) and give the contour plot in Fig. 2 for tan/5 = 50. 
From Fig. 1, we know that the most of the region corresponds to n/2 < a < n. Thus, 
the point moves to the upper right corner in sin 2 26*i 2 and m x plain, a approaches to 
~ 7r from tt/2 and both sin 2 26 l i 2 and mi increase. Since \Vi 3 \(mi/m 3 ) is proportional 
to m 2 sin26>i 2 sina:o/2, its value increases rapidly. This situation is seen from Fig. 2 for 
tan/5 = 50. Thus, we may easily expect the value as large as 0.05. In order to obtain 
I V13 1 , we have to multiply m^/mi which may push its value larger, if m 3 > mi. 

(c) The induced Dirac CP violation phase 5 

The induced Dirac CP violation phase contains /5 which is not fixed in this model. 
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Therefore, in general, the Dirac phase can take any value, until we fix the value of (3 . 
In order to estimate the Dirac phase aside from (3 , we define Si given in Eq. (29), which 
is obtained by excluding /3 . With tan/3 = 50, we show in Fig. 3, values of sin<5i in the 
sin 2 2# 12 and mi plain. The solid line shows a curve on which sin^ takes a fixed value. 
In the right-half domain, the larger value is obtained. As we stated, if a® enters in the 
domain — n < a < 0, the sign of sin^i changes. 

(d) The effective mass of the neutrinoless double beta decay (m v ) 

The effective mass (m u ) is proportional to mi as shown in Eq. (30), it becomes larger 
as mi increases, while decreases if sin 2 2# 12 increases. 

All corresponding figures for tan/3 = 20 are shown in Figs. 4, 5, 6. Except for (m u ), 
the figures are obtained by shifting the larger mi, because the dependence of tan/3 is 
scaled by em\ as we can see from the definition of h in Eq. (24). The effective mass (m u ) 
is almost the same as the case of tan/3 = 50. 

5. SUMMARY AND DISCUSSIONS 

In this paper, we discuss a class of neutrino mass matrix which predicts zero or a 
small value of |Vi3| and found the inequality in Eq. (28). This constraint gives a severe 
restriction for model building of neutrino mass matrix. In particular, the model which 
predicts a larger value of tan 2 6> 12 at M R scale than the experimental value obtained from 
the solar neutrino mixing is excluded. As a result, the bi-maximal mixing scheme at M R 
scale is excluded, if the experimental value tan 2 6 & < 1 is established. 

In this model, | V13 j in Eq. (19) at mz which is induced radiatively may not be small 
as it is shown in Fig. 2, if the neutrino mass mi is of order 0.05 eV. The Dirac phase Si 
in Eq. (29) at m z which is also induced may not be small in general as we see in Fig. 3. 
The effective neutrino mass (m u ) in Eq. (30) is expected to be of order 0.05 eV. All these 
values for |Vi 3 |, Si and (m u ) depend crucially on the mass mi which is assumed to be 
around 0.05 eV. 

The fact that Majorana phases at Mu scale can induce a Dirac phase pushes our 
dream further to consider the possible relation between a Dirac phase which appears 
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FIG. 1: Contour plot of ao in sin 2 29\2 and mi plain to reconstruct the experimental value of 
9@ in the case of tan/3 = 50. We use as experimental values tan 2 6q = 0.34(sin 2 2(9 Q ~ 0.75), 
Am| = 5 x 10~ 5 [eV 2 ], sin 2 26> atm = 1, and Am 2 tm = 3 x 10~ 3 [eV 2 ]. The allowed region is 
between qq = vr and ao = curves. 



in the neutrino oscillations and the Majorana phase which appears in the leptogenesis. 
We believe such scenario does exist and the finding of the missing link will be the most 
wonderful and fruitful project. 

Recently, Antusch et al.||17|| studied the quantum effect for the neutrino mass matrix 
which reproduces the Bi-Maximal mixing at the GUT scale, Mq. They considered the 
quantum effect due to heavy Majorana neutrinos. They considered two cases, (i) the 
standard model (SM) and (ii) the MSSM with tan/3 = 5. In both cases, the quantum 
effect from the lightest heavy Majorana neutrino mass, Mr to mz is very small so that 
it can be neglected. Therefore, they considered a possibility that that tan 2 # so i = 1 at 
Mq reduces to the experimental value tan 2 # so i = 0.34 at Mr by the radiative correction. 
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FIG. 2: Contour plot of |Vi3|(mi/m3) in sin 2 2du and m\ plain for tan/3 = 50. We use same 
values as Fig. 1 for experimental values. Gray curves show the «o values as in Fig. 1. 



They found it possible in a special situation where the Dirac mass matrix is in the form 
of m D = (v /v / 2)diag(l, e', e /2 ) with e' < 1, which in turn means that M R ~ 10 14 GeV. 
This model has two special features: One is that the Dirac mass matrix has the inverse 
mass hierarchy which disagrees with the naive expectation from the GUT scenario. The 
other is that the scale of M R is larger than the ordinary expectation, M R ~ 10 13 GeV in 
order to have larger Yukawa coupling constants, y u , related to the Dirac neutrino mass. 
In our case, we consider the large tan/5 case so that Yukawa coupling constants y u are 
small. Thus, the correction is negligible and our result is valid even at M G . 
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tan(3=50 




FIG. 3: Contour plot of sin<5i and (m„) for tan/3 = 50. We use same values as Fig. 1 for 
experimental values. Solid curves denote sin<5i and dashed curves denote (m„). Gray curves 
show the ao values as in Fig. 1. 



Appendix A: Diagonalization of the neutrino mass ma- 
trix 

We define diag(l, 1, a)0 = OX, 



X = 1 - eO T diag(0, 0,1)0 



= 1 



(22 2 ^ 

S 12 S 23 s 12Cl2S 2 3 Sl2S23 c 23 

2 2 2 

S12C12S23 C 12 S 23 c 12 s 23 c 23 



(A.1) 



y-Sl2S23C23 C12S23C23 C 23 J 

where e = 1 — a is a small positive quantity and its value is given in Eq. (8), then we 
consider the mass matrix transformed by O as 



m u = T m u (m z )0 = XD V X A 



(A.2) 
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tan(3=20 

sin 2 26 12 

. 7 
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0.02 0.04 0.06 0.08 0.1 0.12 

mJeV] 

FIG. 4: Contour plot of ao in sin 2 26? 12 and mi plain to reconstruct the experimental value of 
6 & in the case of tan/3 = 20. We use same values as Fig. 1 for experimental values. The allowed 
region is between ao = ir and ao = lines. 

Now, we diagonalize m„. 

In order to diagonalize this matrix directly, we consider the Hermite matrix m\,m v 

m\m v ~ (1 - 4e S 2 2 s 2 3 )M 2 + Y , (A.3) 

where elements of Y are given up to the 1st order of e as 



Y u 


= , Y 22 


— A21 , y 33 — A31 , 


Y 12 


— r 21 — " 


-2esi 2 ci 2 sl 3 (M 1 2 + M| + 2M 1 M 2 e ia °)M 1 2 


Y 13 


V* 

— *31 ~~ 


-«i2S23C23(M 2 + M 3 2 + 2M 1 M 3 e ,ft ) , 


Y23 


— 2 32 — 


-ec 12 s 23 c 23 (M 2 + M 2 + 2M 2 M 3 e l(ft " o) ) 



(A.4) 
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FIG. 5: Contour plot of |Vi3|(mi/m3) in sin 2 2#i 2 and mi plain for tan/3 = 20. We use same 
values as Fig. 1 for experimental values. Gray curves show the «o values as in Fig. 4. 



A 21 = A 21 - 4es 2 3 (4M 1 2 - s 2 2 M 2 ) , 

A 31 = A 31 - 4e(c 2 3 M 3 2 - s 2 12 s 2 23 M?) . (A.5) 

Since mj,m„ is an Hermite matrix, it is diagonalized by the unitary transformation as 
V^mlm u V ~ (1 - 4es 2 2 s 23 )M 2 + V^YV. The diagonalization of the matrix Y can be 
achieved by using the see-saw technique, because | ^33 1 is much larger than all other terms. 
That is, by using the unitary matrix 





/ 


1 





Y 13 \ 








Y33 









1 


Y 23 
Y33 




V 


'l3 
Y33 


^23 
133 


V 
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FIG. 6: Contour plot of sin^i and {m u ) for tan/3 = 20. We use same values as Fig. 1 for 
experimental values. Solid curves denote sin<5i and dashed curves denote (m„). Gray curves 
show the ao values as in Fig. 4. 



Y is block diagonalized in a good accuracy as 

/ _ 2|y- 13 | 2 v 2) L ,)-, 



V^YV 3 



Y 33 

_ 2y i*3 y 23 

r 12 Y 33 



Y 



12 



^22 — 



Y 33 

2|y 2 3l 2 

Y 33 



/ Y 12 ^ 



Y* 2 Y 22 
Y 33 



^33 + 






\y 13 \ 2 +\y 23 \ 2 

Y 33 



\ 



(A.7) 



where in the last equation, we neglected the see-saw induced terms because \Yi 3 / Y 33 \ and 
1^23/^33! are much smaller than 1 and Yi 3 and Y 23 are same order of Y i2 and Y 22 . 

In the following, we use Mi = M 2 for terms proportional to e to simplify the expression. 
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The matrix in Eq. (A7) is diagonalized by 





( 


c 


_ se ia /2 


o\ 


v l2 = 


se 


-iao/2 


C 







V 








V 



with c = cos 9 and s = sin 9 as 

^12(^3^3)^12 = diag(Ai, A 2 , A 3 ) , 



where 



tan 29 = 



4eM 2 sin2# 12 s 2 3 cos ^ 



V21 



and 



Ai 



^21 



1 



, A 2 



^21 



1 



, A 3 



2 V cos 20/ - 2 V cos 29 
Neutrino masses at m z are obtained by m 2 = (1 — 4es 2 2 s 2 ,3)M 1 2 + A; and 

A 2 i 



31 



A ^atm 



m 9 — mi = 



1 cos 20 ' 

2 _ _2 2 



m 2 — m 3 — m 1 ~ A31 . 



We find that the angle # is expressed by 



sin 2^ 



4eM? sin 20i 2 s| 3 cos f 



Am| 



By taking into account of 



(OV 3 V 12 ) T m u (m z )OV 3 V 12 ~ diag(mi, m 2 e ia °, m 3 e <A) ) , 



(A.8) 



(A.9) 



(A.10) 



(A.11) 



(A.12) 



(A.13) 



(A. 14) 



with 777^ > 0, we find the mixing matrix V which satisfies V T m v {mz)V 
diag(mi, 777,2, m 3 ) is V = OV 3 l / i 2 diag(l, e~ ia0//2 , e~ l/3 °/ 2 ) is given by 

'12' ' " '" V>"> 12' " '" ! 2-"" ' — I ' 13|< "'^ ' ' '' ' ' 

V ~ 



c 12 c- s 12 se~ tao ' 2 -(s 12 c + c 12 se ia °< 2 ) -\V 13 \e tp 
c 23 (si 2 c + c 12 se- ia °/ 2 ) c 23 (ci 2 c - s 12 se ia °/ 2 ) -s 23 
\s 23 (s 12 c + c 12 se- ia °/ 2 ) s 23 (c 12 c-s 12 se ia °/ 2 ) c 23 ) 



1 

e~ ia °/ 2 







-iPo/2 



) 

(A.15) 
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